LyapunovSolve
LyapunovSolve[a,c] finds a solution x of the matrix Lyapunov equation a.x+x.ac.
LyapunovSolve[a,b,c] solves a.x+x.bc.
LyapunovSolve[{a,d},c] solves a.x.d+d.x.ac.
LyapunovSolve[{a,d},{b,e},c] solves a.x.e+d.x.bc.
Please visit the official Wolfram Language Reference for more details and examples on core symbols.